Long wavelength structural anomalies in jammed systems 
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The structural properties of static, jammed packings of monodisperse spheres in the vicinity of the 
jamming transition are investigated using large-scale computer simulations. At small wavenumber 
k, we argue that the anomalous behavior in the static structure factor, S{k) ~ k, is consequential of 
an excess of low-frequency, collective excitations seen in the vibrational spectrum. This anomalous 
feature becomes more pronounced closest to the jamming transition, such that S{0) ^ at the 
transition point. We introduce an appropriate dispersion relation that accounts for these phenomena 
that leads us to relate these structural features to characteristic length scales associated with the 
low-frequency vibrational modes of these systems. When the particles are frictional, this anomalous 
behavior is suppressed providing yet more evidence that jamming transitions of frictional spheres 
lie at lower packing fractions that that for frictionless spheres. These results suggest that the 
mechanical properties of jammed and glassy media may therefore be inferred from measurements of 
both the static and dynamical structure factors. 

PACS numbers: 61.43.Fs64.70.Pf81.05.Rm 



I. INTRODUCTION 

The emergence of similarities between the properties 
of molecular glasses, dense colloidal suspensions, foams, 
and granular materials, have led to the notion of jamming 
[H - the transition between solid-like and fluid-like phases 
in disordered systems - as a manner through which one 
can gain a deeper understanding of the traditional liquid- 
glass transition and the fascinating, complex phenomena 
observed in amorphous materials in general ■ Although 
there has recently been some works highlighting the dif- 
ferences between the onset of rigidity in jammed matter 
and glassiness 0, 0] , the emphasis in this work is on the 
commonalities between the two 0, @, 0, B- Here, we aim 
to provide a heuristic physical picture that accounts for 
specific, long wavelength, structural features that emerge 
in the jammed state and relate these features to their dy- 
namical properties. 

Donev et al. found that in the hard sphere jamming 
transition, the structure factor, S{k), exhibits a linear 
dependence on wavenumber k = |k|. 
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This behavior of S{k) suggests that the total correlation 
function h{r), decays as l?""^!, at large separations r, as 
deduced from the asymptotic estimates of Fourier trans- 
forms [13, M, El- Similarly it suggests a long range 
behavior for the direct correlation function c{r). This 
is indeed in contrast to standard liquid state theory for 
liquids whose constituents interact via a finite range po- 
tential - hence c(r) is short ranged - which predicts, 
S{k) (X k^. This anomalous low-fc behavior can be in- 
terpreted as being indicative of the suppression of long 
wavelength density fiuctuations due to hyperuniformity 
p^ . Here, we propose an alternative interpretation re- 
lated to large length-scale collective dynamics. 

This paper is arranged as follows. We provide a brief 



overview of the Molecular Dynamics (MD) simulations 
used here to generate liquid and jammed states. We 
then review previous results from studies of the jamming 
transition pertinent to the discussion here. This is fol- 
lowed by a discussion on the relevant concepts from liquid 
state theories that indicate that our results and those of 
Ref. 9 for frictionless particles are indeed rather unusual. 
We then present our results for the static structure fac- 
tor, S{k), at small values of the wavenumber fc, in our 
jammed, model glassy system. We then put forward a 
conjecture that relates the asymptotic behavior of S{k) 
to an excess of vibrational modes relative to the Debye 
model. We end with results from ongoing work on fric- 
tional systems and conclusions. 



II. SIMULATION MODEL 

The computer simulations performed here are for 
monodisperse, soft spheres of diameter d and mass m, 
interacting through a finite range, purely repulsive, one- 
sided, harmonic potential. 
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where r = jr^ — rj\ is the center-to-center separation be- 
tween particles i and j located at r^. The strength of 
the interaction is parameterized by e, which is set to 
unity in this study. Most of the results presented be- 
low are for frictionless particles, with friction coefficient 
fi — 0. We also present preliminary results for frictional 
packings using a static friction model [3] to compare 
between frictionless and frictional systems. In the fric- 
tional packings the particle friction coefficient was varied, 
0.01 < /i < 1.0. We simulated systems ranging in size, 
1024 < N < 256000 particles, in cubic simulation cells of 
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size lOd ^ L 50d, with periodic boundary conditions, 
over a range of packing fractions, = Nird^ /L^. 

The jamming protocol implemented here to generate 
zero-temperature jammed packings is similar to other 
soft-sphere protocols [IB, . Initially, starting from a 
collection of spheres randomly placed in the simulation 
cell at low packing fraction 0^ — 0.30, we compressed the 
system to a specified over-compressed state, = 0.74, 
minimizing the energy of the system in a steepest de- 
scent manner. At this value of = 0.74, all the particles 
experience overlaps with several other particles - their 
contact neighbors. The packings are mechanically sta- 
ble and disordered with no signs of long range order. To 
generate packings at < 0.74, we then incrementally 
reduced in steps of, 10~^ < 6(l> < 10"^, minimizing 
the energy after each step. This allowed us to accurately 
determine the location of the jamming transition where 
the system unjams at a packing fraction 0c (see below), 
for each independent realization, down to an accuracy of 
10^^ in (j) for N < 10000 and an accuracy of 10^^ for 
N > 10000. For the largest system, N = 256000, we 
generated four independent configurations at each value 
of 0. 

We also studied equilibrated liquids of frictionless 
spheres at a dimensionless temperatures T = 0.01. This 
value of T was chosen to compare and contrast our results 
between liquid and jammed states. 

Our main tool of analysis is the static structure factor 
S{k), which we obtain in the standard way as a direct 
Fourier transform of the particle positions. 
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Here, k is the reciprocal lattice vector for the periodic 
simulation cell, which is restricted to fc > 27r/L. Al- 
though the data is noisy at the lowest wavenumbers, our 
largest packings {N = 256000) display the anomalous 
low-fc behavior, as described below, well above the noise 
level. Throughout this study, d = m — e = 1, and all 
quantities have been appropriately non-dimcnsionalizcd. 



III. REVIEW 

A. Jamming Transition 

Frictionless, purely repulsive, soft sphere system have 
been shown to undergo a zero-temperature transition be- 
tween jammed and unjammed phases as a function of 
density or packing fraction [B, [3, [l^, [13, [H, [l^. In 
the infinite-size limit, this transition occurs at a packing 
fraction coincident with the value often quoted for ran- 
dom close packing, (j)rcp = 0.64 Above the jamming 
threshold > 0c, the jammed state is mechanically sta- 
ble to perturbations with non-zero bulk and shear mod- 
uli, whereas below the transition < 0c, it costs no 



energy to disturb the system [1^. The relevant param- 
eter here is not the absolute value of the packing frac- 
tion, rather the distance to the jamming transition de- 
fined through, A0 = — 0c. Thus, as a jammed state is 
brought closer to the transition point, A0 — > 0, it gradu- 
ally loses its mechanical rigidity and becomes increasingly 
soft. Intriguingly, this loss of mechanical stability as the 
jamming transition is approached can be related to di- 
verging length scales Q that characterize the extent of 
soft regions that determine the macroscopic behavior of 
system j23|- Experiments [2lj have also identified grow- 
ing length scales in the vicinity of the jamming transition 
thus promoting the idea that the jamming transition can 
be considered in the context of critical phenomena. 

In this paper, we connect the long wavelength 
structural features observed in S{k) to correlation 
lengths characterizing the typical length scale of collec- 
tive, low-frequency vibrational modes in jammed, zero- 
temperature, disordered packings Jj. In Ref. |7l it was 
shown that the approach of the jamming transition in 
a soft sphere packing is accompanied by a dramatic in- 
crease in the number of low-frequency vibrational modes 
over the expected Debye behavior. In traditional glasses, 
these excess low-frequency modes in the vibrational den- 
sity of states, 'D{oj), are often referred to as the boson 
peak "22*1 in reference to the peak observed when plot- 
ting T>{lu)/u!'^. For convenience, we also employ this lan- 
guage here. A detailed study of D{uj) for jammed sphere- 
packings and the appearance of the so-called boson peak 
can be found in Q . 

In Fig. [T] we compare the vibrational density of states 
for dense, soft sphere liquids (Fig. [TJa)) and amorphous 
jammed solids, identifying the location of the boson peak, 
lob, for our jammed system (Fig. IHb)). The location of 
the boson peak tends to zero at the jamming transition 
point, i.e. ujb 0, as A0 — > Q. The two values of 
A0 = 1 X 10"'' and 1 x 10^^ correspond to actual packing 
fraction values of = 0.6405 and 0.74 respectively, for 
sample configurations used to generate this data. 

Associated with this excess in the vibrational density 
of states are two diverging length scales: the longitu- 
dinal correlation length, ^i, characterizing the scale of 
collective excitations of longitudinal modes contributing 
to the boson peak, while ^t, characterizes transverse ex- 
citations. These correlation lengths scale with the boson 
peak position as 0, H^l , 

ujb oc (L: longitudinal modes) (4) 

cjb oc ^j^^ (T: transverse modes). (5) 

Thus, when the Debye contribution, lod, is included the 
corresponding dispersion relations become, 



ujL{k) - ujd,l = Pk 
LJT{k) ~ lod.t — ak^ . 



(6) 
(7) 



The preceeding relations are only approximate, valid for 
low frequencies. They both contribute to the boson peak: 
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FIG. 1: The vibrational density of states, "Diuj), for two sys- 
tems of A'^ = 1024 soft, frictionless spheres, characterized by 
the distance = (j) — cj)c, from their zero-temperature, jam- 
ming transition packing fraction. Solid line: A.<f) = 1 x 10~^. 
Dashed line: = 1 x 10~^. (a) Equilibrated liquids 

at T = 0.01, well above their respective freezing or ap- 
parent glass transition temperatures. Data obtained from 
Fourier transforms of velocity autocorrelation functions, and 
has been rescaled such that the u) = Q intercepts coincide. 
A Lorentzian function (dotted) corresponding to Langevin 
diffusion [23| is also shown, (b) Jammed packings at zero 
temperature. The location of the boson peak is identified as 
UJB, for A(j> = 1 X 10"''. In (b), the Debye result (dotted), 
corresponding to the system A(j) = 0.1, was generated using 
values for the bulk and shear moduli 6], and is shown for 
comparison. 



The longitudinal term by modifying the slope in the De- 
bye relation, whereas the transverse contribution con- 
tains an anharmonic term. We discuss the relevance of 
these results further below. 



B. Analyticity of S{k) 

In order to avoid any misunderstanding we define the 
concepts of regular and singular in the following sense 
[To| . Assuming A; to be a complex variable, for any three 
dimensional system interacting with tempered pairwise, 
additive potentials, namely when 
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as r 
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and therefore integrable, |/ dr(l){r)\ < oo. Then, the 
structure factor, in the long wavelength limit reads [l^ . 



Sik) : Aplfcl" -HF(fc2) 



as 
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where p = y is the number density, A is a. constant that 
depends on the thermodynamic properties of the equilib- 
rium system, and F is an analytic function of fc, regular 
in the neighborhood of fe = 0. Thus, the second term 
on the right hand side of Eq. [9] is the regular contribu- 
tion to S{k), whereas the first term - originating from 
the potential of interaction - is the singular term. 

For instance, when r; = 3, as in the attractive part of 
the Lennard- Jones potential, 



S{k)^Ap\k\^ + F{k'^) 



(10) 



a result originally derived by Enderby, Gaskill, and 
March Thus, a linear behavior of the structure fac- 
tor at small values of k corresponds to 77 = 1. Physi- 
cally, this would correspond to a charge-dipole interac- 
tion, which Chan et al. [25*1 have shown can only be 
present if long-ranged dipole-dipole interactions are also 
included. Thus, in the present work, the linear behavior 
of 5'(fc), at small fc, does not come from the potential of 
interaction, as the finite range harmonic potential used 
here will only give the regular contribution to S{k). How- 
ever, this does suggest that the anomalous low-fc behavior 
originates in long-ranged correlations associated with the 
system. 



IV. RESULTS 

A. Frictionless Spheres 

We initially generated a number of over-compressed, 
4> > 4>c, zero-temperature packings at various distances, 
A(f) = (f) — (f)c, from the jamming transition point. We also 
point out that the main distinction between this work 
and that of Ref. 9, is that we study soft-spheres and 
approach the jamming transition from above, (j) — > c/)'^ , 
whereas, hard spheres necessarily approach the transi- 
tion from below. In Fig. [2] we show S{k) for two soft- 
sphere systems at T = 0.01, well above their respective 
apparent glass transition temperatures Tg. (The corre- 
sponding density of states are shown in Fig. m^a).) For 
A(j) = 10-\ Tg < 0.001, while for smaller A(/>, T„ lies 
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FIG. 2; Two iV = 10000 soft-sphere systems in the Uquid 
state, T = 0.01, for A<^ = 1 x 10"^ and 1 x 10"* (see legend 
in panel (b)). This temperature is well above the respective 
glass transition temperatures for the two systems, T » Tg. 
(a) Static structure factor S{k), as a function of wavenumber 
k. The inset is a zoom of the region near the origin, where 
the arrows indicate the non-zero intercepts as fc 0. (b) 
Radial distribution function, g{r) shows a prominent nearest 
neighbor peak at r ^ 1 and correlations that die off rapidly. 



below this. Qualitatively the curves are similar. We 
find the usual primary peak corresponding to nearest 
neighbors, and, as expected from standard liquid state 
theory [H, HI], at lower fc, S{k) oc fc^, with a non- 
zero intercept at = 0, which corresponds to the fi- 
nite compressibility sum rule, Fig. [2] inset. Thus, for 
equilibrated liquids the structure factor and the long- 
wavelength limit are quite insensitive to the location of 
the zero-temperature jamming transition, at these values 
of the number density, p = N/L^ = Gcj^/ncf = 1.22, 1.44 
for = 0.6405 (A0 = 1 x 10-*), 0.74 (A0 = 1 x IQ-^) 
respectively. 

Turning our attention to the jammed phases at zero 
temperature [26| , in Fig. [3] we show the radial distri- 
bution function g{r) and the low-A; region of S{k), for 
N — 256000 soft-spheres systems at three values of 
A(j) = 1 X 10-\ 1 X 10-2, and 3 x lO'^. Note that in 



Fig. [Slja) we use log-log scales to clearly demonstrate the 
linear region in S{k). Far from the jamming transition, 
A(p = 1 X 10^^, S{k) plateaus near fc « 1 and tends to 
constant as fc ^ 0. For the systems closer to the jam- 
ming transition. At/) = 1 x 10"^ and Ac/) = 3 x 10^"^, 
there is a qualitative change in the low-fc behavior of 
S{k). We find a linear region, S{k) ^ fc, extending over 
almost an order of magnitude at low-A:. We also point 
out that on closer inspection the linear region extends to 
lower k for the system closest to the jamming transition 
point (A0 = 3 X IQ-'^). At the smallest k attainable, 
S{k), flattens out again because the system is not ex- 
actly at A(j) = 0. The radial distribution function shown 
in Fig. ^h) exhibits typical features characterisitc of a 
glassy phase; namely a split second peak an additional 
shoulder on the third peak. 

For completeness, we examine the influence of system 
size on the results presented here. In Fig. 21 we show S{k) 
at A(f) « 3 X 10"^ for three different systems sizes, N = 
1000, 10000, 256000, corresponding to L « 10, 20, 60. 
The main panel of Fig. 2] shows that the gross properties 
of S{k) do not depend on system size. Oscillations in 
S{k) persist out to the largest k and is a consequence 
of the diverging nearest neighbor peak in g{r) Q. The 
inset to Fig. 2] indicates that the linear portion of S{k) 
at small k becomes resolvable for N > 10000. 

Thus the flrst aim of this work shows that, similarly to 
the flndings in Ref. for hard spheres just below the jam- 
ming transition point, for our soft sphere packings above 
(/>c, S{k) exhibits a linear behavior at small values of mo- 
mentum transfer. Hence, this unusual behavior not only 
pertains to hard spheres, but also to soft spheres, both 
above and below the jamming transition point, and pos- 
sibly, to flnite range repulsive potentials in general. We 
will attempt below to produce what we believe is a rea- 
sonable explanation for the origin of this behavior which, 
as stated at the beginning, does appeal to a conjecture. 

Since, in the long wavelength limit, the linear behavior 
of S{k) cannot arise from the singular contribution of the 
potential used in this work [l^, [3| , it can only be due to 
the collective excitations present in the jammed or glassy 
state [131 • The collective excitations we have in mind are 
the same that are responsible for the boson peak. The 
range of fc, over which S{k) ~ fc, characterizes the length 
scale over which these excitations may be considered col- 
lective. As one approaches the jamming transition, this 
length scale diverges, a result consistent with work on 
the density of states of jammed packings 

Recent work by Chumakov et al. [28*1 does indeed show 
that the excitations leading to the boson peak are pre- 
dominantly collective; in agreement with inelastic neu- 
tron scattering experiments Although there have 
been a number of theories put forward [30., 31, 32], there 
is, at present, no agreed explanation as to the origins of 
the boson peak. For the purposes of this work this is 
not necessary, except insofar as there is agreement that 
these excitations are predominantly collective, and that 
the vibrational modes are likely to be kinetically driven. 
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FIG. 3: Zero-temperature, jammed packings containing A'^ = 
256000 purely-repulsive, frictionless, soft-spheres at three dif- 
ferent values of A0 = 1 x 10"^ (dotted line), 1 x 10"^ 
(dashed), and 3 x 10"'^ (solid), (a) Static structure factor 
S{k). Far from the jamming transition, Acji = 1 x 10~^, S{k) 
plateaus nears « 1. Closer to jamming, A(/> = 1 X 10"^ and 
Atj) = 3 X 10"'^, ^(A:) exhibits approximately linear depen- 
dence on k over almost an order of magnitude in k, extending 
down to low-fc. A linear curve on this log-log plot is shown 
for comparison, (b) Radial distribution function. The jam- 
ming transition is characterized by a diverging nearest neigh- 
bor peak at r ~ 1, a clear splitting of the second peak Q], 
and oscillations that persist out to larger r than for the liquid 
state. 



The arguments that follow are somewhat oversimpUfied, 
but we believe they are along the correct lines. 

The following comments are in order: (i) In a system 
that undergoes the liquid-glass transition, the long wave- 
length limit of the static structure factor, S{k — > 0), is 
not the thermodynamic compressibility kt, but a differ- 
ent value, sa y k , which can actually be extracted from 
experiments [33|. (ii) Independently of the model or ap- 
proximation used, the density fluctuations in the glas s 
phase are kinetically, not thermodynamically, driven [33| . 
(iii) In spite of their different origins, the kinetic temper- 
ature, say Tdyn, defined for glasses, matches the Edwards 
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FIG. 4; S{k) at A(l> ^ 3 x 10~^ for different systems sizes, 
iV = 1000 (line), iV = 10000 (o), and TV = 256000 (+). Oscil- 
lations persist out to large k, reflecting the dominant nearest 
neighbor peak in g{r) shown in Fig. [3] The inset indicates 
that the linear behavior at low-fc becomes better resolved with 
increasing system size. 



temperature, T^dw, used in granular matter. They hap- 
pen to coincide within mean field theories when the glass 
is in contact with an almost zero temperature heat bath, 
and the athermal grains are jammed [34 1. 

The connection between the dispersion relations and 
the structure factor, follows from the second moment of 
the dynamical structure factor, S{k,uj) [35| . 
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dujLU S{k,uj) — VQk 



(11) 



where z/q — vq{T) (for a system in thermodynamic equi- 
librium Vq — -^^j where ks is Boltzmann's constant). 
The static structure factor is the zero moment of S{k,uj), 



S{k) = / dujS{k,uj). 



(12) 



We now put forward the following conjecture: assume 
that there is only one very well defined collective mode 
with dispersion relation ojB{k), the boson peak. This 
gives us a relation between the asymptotic long wave- 
length behavior of S{k) and the boson peak. 



whence, 



S{k,u}) ^ S{k)6{Lo - ujB{k)) 



S(k) = 4 



2^1 



(fc) 



(13) 



(14) 



This collective mode is associated with vibrational exci- 
tations with a wavelength of the order of the correlation 
length of the jammed state. We expand the dispersion 
relation at these small values of momentum transfer 12711 



ujB{k) = ck + ak^ + 



(15) 
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In Eq. [151 c is the speed of sound, and the second term 
on the rhs denotes departures from the usual Debye be- 
havior, such that both c and a are independent of k. 
Replacing Eq. [15] into [TH in the limit of fc ^ 0, we find, 

Thus, on comparing Eqs. [TSj and [TBj to Eqs. [6] and [3 it 

transpires that it is the transverse modes that contribute 
to the linear behavior of S{k). On the other hand, the 
longitudinal components, Eq. [B] only contribute a con- 
stant term to S{k), modifying the slope of the Debye 
relation. 

We illustrate the differences in the dispersion behavior 
at two extreme values of A(p — 10~^ and 10~^ in Fig. [5] 
These results were obtained from diagonalization of the 
dynamical matrix for systems with N — 10000 and then 
locating the peaks in the transverse components of the 
Fourier transforms of the eigenmodes - transverse mode 
structure factors [s^ [37j - and then averaging over a 
small range of frequencies. The data shown in Fig. [5| dis- 
tinguishes between the regular, linear dispersion relation 
that dominates the dispersion behavior far from the jam- 
ming transition at At/) = 10"^, while the quadratic con- 
tribution to the dispersion behavior is significant closer 
to the jamming transition, A0 = 10^^. 

1 I . . . . . . . 1 




FIG. 5: Transverse dispersion behavior for A(j) = 10~^(o) 
and 10~^(n). The solid lines correspond to quadratic fits to 
the data as in Eq. \16\ Data obtained from the low-frequency 
portion of the transverse structure factors of the vibrational 
modes for N = 10000. 



From the analysis we obtain the transverse speed of 
sound, Cf, as a fitting parameter. At these two com- 
pressions: ct « 0.26(A(/) = 10-1) and a ~ O.O23(A0 = 
lO"®). The actual values calculated from the bulk and 
shear moduli data Q give: Cj « O.28(A0 — 10^^) and 
Cf ~ 0.018(A(/) = IQ-^). Therefore, we find reasonable 
agreement within this approximation. If we use only a 
linear dispersion relation we find that the values of ct 



computed here and the actual values [^ vary by as much 
as an order of magnitude. This can be realized by the ob- 
servation that a linear fit to data for A(p ~ 10^^ does not 
pass through the origin as required in the hydrodynamic 
limit. 

The linear feature of S{k) in the small-fc regime ap- 
pears here as a consequence of assuming an "excess" rel- 
ative to the Debye model for the dispersion relation, and 
only one collective mode. Although the "excess" model 
used here appears as an approximation appropriate for 
the low-fc regime, we have reconciled this with existing 
results on the emergence of characteristic length scales 
associated with these collective excitations. More gener- 
ally, however, it does suggest that the boson peak and 
the linear behavior of S'(fc), as k — > 0, are two sides of 
the same coin. Therefore, this anomalous suppression of 
long- wavelength density fluctuations emerges as a conse- 
quence of large length-scale correlated dynamics in the 
low-frequency modes of the jammed solid. Furthermore, 
although the boson peak has been traditionally associ- 
ated with the glassy phase, recent low frequency Raman 
spectroscopic studies of glassy, supercooled, and molten 
silica reveal that the boson peak persists into the liquid 
phase [sll. This also appears to be the case with other 
network systems (see references in Ref. Issh which, in the 
molten state, show a distinctive prepeak at small k in 
S{k). This prepeak is indicative of intermediate range 
order in those melts, representing a characteristic length 
in those melts that measures the correlations between 
the centers of "clusters" present in the liquid state. This 
intermediate-range clustering may in fact promote longer 
wavelength correlated dynamics in the liquid state. It 
may therefore be interesting to investigate experimen- 
tally, and by means of MD simulations, whether the low 
fc-behavior of S{k) in those systems is also linear, both 
in their glassy and liquid phases. 



B. Prictional Packings 

Studies on jammed packings of frictional particles are 
less well-developed. The picture that is emerging is that 
frictional packings undergo a similar zero-temperature 
jamming transition at packing fractions which now be- 
come friction dependent [H, [40, [HI, . In the limit of 
high friction coefficient the frictional jamming transition 
coincides with the value associated with random loose 
packing, 0^;^ « 0.55 [H, jH, Hi] . Thus the relevant pa- 
rameter that measures the distance to the jamming tran- 
sition now becomes a friction-dependent quantity, A(/)(/i) 

We test these concepts through S{k) for packings over 
a range of friction coefficients at the same = 0.64, 
which translates to different values of A0(/Lt). Our pre- 
liminary results in Fig. [6| are data for S{k) for jammed 
packings with < n < 1. At this fixed value of 
(j) = 0.64, A4>{fi) increases with increasing friction co- 
efficient. Therefore, we expect that the anomalous linear 
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FIG. 6: Log-log plot of the static structure factor S{k), 
focusing on the low-fc region. All packings at the same 
(j) = 0.64, contain A'^ = 256000 purely-repulsive, monodis- 
perse, soft-spheres, with different friction coefficients /i = 
(solid - line), 0.01 (dash), O.l(dot), 1. (dot - dash). 



behavior in the low-fc region of S{k) should become less 
prominent with increasing fi. Indeed this trend is ob- 
served in Fig. [6l thus providing yet further evidence that 
the jamming transition point shifts to lower </> with in- 
creasing friction. To recover the linear behavior at low-A: 
in frictional packings we would therefore need to study 
mechanically stable packings at lower 0. This work in 
ongoing. 

V. CONCLUSIONS 

In conclusion, we have provided a physical, albeit 
naive, explanation for the observation of the apparent 
suppression of long wavelength density fluctuations in 
jammed model glassy materials. We submit that it is the 
presence of low-frequency, collective excitations, mainly 
of transverse character, contributing to the excess of low- 
frequency modes, that are responsible for the linear be- 
havior at low-fc in S{k). The relevant length scale here is 
the transversal correlation length that contributes to the 
diverging boson peak in the jammed phase. This low-fc 
feature is most pronounced in the model system studied 
here for soft particles at packing fractions above the zero- 
temperature jamming transition matching the behavior 



on the other side of the jamming transition for hard 
spheres. This connection between the long- wavelength 
behavior on cither side of the transition reinforces the 
view that the jamming transition is critical in nature. 
Although we also point out that the jamming transition 
is very different from a thermodynamic liquid-gas critical 
point in that we do not observe a divergence in ^(fc — 0), 
but rather it becomes zero. 

Moreover, we expect that as the transition is ap- 
proached, 

i^o(T) -> ^ S'(O) ^ as A(t)-^0. (17) 
The dependence of t^oiT) on a generalized temperature- 
like quantity, T, in jammed packing s is consistent with 
the concept of the angoricity |46l . |47[ . This plays the role 
of temperature in packings of elastic particles and goes to 
zero when the particles fall out of contact which occurs at 
the jamming transition. Future studies may also provide 
a way to understand the connection between different 
temperature definitions relevant to the study of thermal 
and jammed systems. Our preliminary data for frictional 
packings also provide further evidence that the location 
of the jamming transition occurs at lower packing frac- 
tions with increasing friction coefficient. The underlying 
nature of the low-frequency modes in frictional materials 
has yet to be investigated in three dimensional systems 
and forms part of ongoing work. When friction is present 
the modes will contain not only the translational char- 
acter as seen in frictionless systems, but also rotational 
character due to the additional degrees of freedom. 

This anomalous low-fc behavior is also likely to be 
present in other, finite range, model repulsive systems. 
These features may be detected, using the appropriate 
spectroscopy, in (hard-sphere) colloidal glasses, and gran- 
ular packings, which can be prepared close to their re- 
spective jamming transitions. 
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